Axino dark matter from Q-balls in Affleck-Dine baryogenesis 
and the flh — floM coincidence problem 
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We show that the 0.^ — Q.dm coincidence can naturally be explained in a framework where axino 
is cold dark matter which is predominantly produced in nonthermal processes involving decays of 
Q-balls formed in Affleck-Dine baryogenesis. In this approach, the similarity of Q,}, and ^dm is a 
direct consequence of the (sub-)GeV scale of the mass of the axino, while the reheating temperature 
Tr must be low, some 10^ GeV, or less. 



o ■ 

o 

o 

(N 

Oh 

<\ 

^ : 

(N . 
(N • 

> : 

cn . 
o 

00 

o 

o : 

i!^ ■ 
Oh: 

Oh 



X 



1. Introduction. 

The origin of nonbaryonic cold dark matter (DM) and 
of baryon asymmetry in the Universe are among the 
longest lasting puzzles in cosmology as well as in par- 
ticle physics today. In particular, the question of why 
the observed values of baryon density f2f, and of dark 
matter VLdm are so close to each other, Q.uM/^b — 
5.65 ± 0.58 [1], remains a mystery. 

A standard paradigm is that the nonbaryonic cold dark 
matter is made up of some weakly interacting massive 
particle (WIMP) which freezes out of thermal equilib- 
rium in the early Universe. Perhaps the most popular 
WIMP candidate is the lightest neutralino of the Min- 
imal Supersymmetric Standard Model (MSSM) as the 
lightest supersymmetric particle (LSP). It remains sta- 
ble due to the conservation of i?-parity. This economical 
scenario does not, however, explain the proximity of Qb 
and fluM- 

The same is generally true for conventional mecha- 
nisms of baryogenesis or leptogenesis. This may indicate 
that the observed baryon-to-DM density ratio is just a 
pure accident, or else a result of some underlying, and 
as yet unknown, more fundamental theory. An alterna- 
tive approach is to try to identify a physical mechanism 
which would simultaneously produce both baryon asym- 
metry and DM in the proportions consistent with obser- 
vations. It is clear that this basically necessitates aban- 
doning standard paradigms for producing both types of 
species in the Universe. This may be one important les- 
son to learn from these considerations. 

A number of attempts at explaining baryon-to-DM ra- 
tio have been suggested in the literature. For instance, 
recently a right handed sneutrino [2] and a sneutrino con- 
densate as an AD field [3] have been proposed. 

A few years ago, Enqvist and McDonald (EMD) pro- 
posed [4, 5] an attractive solution based on a variant 
of Affleck-Dine (AD) baryogenesis [6]. In that scenario, 
an AD condensate forms during inflation and develops a 
large vacuum expectation value (VEV) along a D-flat di- 
rection in the MSSM. Z?-flat directions are configurations 
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of scalar fields for which the Z?-part of the potential van- 
ishes. Prevalent in theories with many scalar fields like 
the MSSM, they are of much interest to cosmology [7]. 

In the standard AD scenario, after the end of inflation, 
the scalar field condensate slowly rolls towards the origin 
and, after a few dozens of coherent oscillations, produces 
a nonzero baryon number in presence of baryon number- 
violating couplings of the fields making up the fiat di- 
rection. Originally, Kusenko and Shaposhnikov argued 
that the AD condensate can instead fragment into non- 
topological solitons called Q-balls [8]. If their baryonic 
charge is large enough, as in models with gauge medi- 
ated SUSY breaking, Q-balls remain effectively stable 
until today, and contribute to the DM density, despite 
severe astrophysical constraints [9]. On the other hand, 
EMD demonstrated that, under nontrivial but natural 
conditions (that we summarize below), in a large class 
of supergravity (SUGRA) models with gravity mediated 
SUSY breaking (GRMSB) Q-balls subsequently decay 
into baryonic matter and neutralino WIMPs assumed to 
be the LSP [4]. 

In the EMD scenario, the baryon-to-DM ratio can eas- 
ily be estimated to be in the right ballpark, as we shall 
see below. This otherwise attractive framework suffers, 
however, from a serious problem: neutralino production 
in Q-ball decays is in fact too efficient, and density 17^^, 
can only agree with observations for low neutralino mass 
~ 1 GeV, well below LEP limits [10]. Moreover, this 
puts into a potential jeopardy the AD mechanism in a 
large class of GRMSB supergravity models. 

In this Letter, we suggest a way out from the above 
problems of the EMD scenario which at the same time 
preserves its successful features, in particular, an expla- 
nation of the r^jj/r^DM ratio. We propose that the DM is 
not made up of the neutralino but instead of an axino, a 
superpartner of the axion. The axino is a neutral Majo- 
rana, chiral fermion. It arises in SUSY models incorpo- 
rating a Peccei-Quinn solution to the strong CP problem 
in QCD. Unlike for the neutralino or gravitino, its mass is 
strongly model dependent and can be much smaller than 
the (gravity mediated) SUSY breaking scale [11, 12, 13]. 
Similarly to the axion, its interactions are suppressed by 
the PQ scale fa — 10^^ GeV, well below the sensitivity of 
LEP. The axino has a number of properties which make 
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it a promising candidate for cold dark matter [14, 15]. 
Earlier papers considered warm axino relics [13, 16]. As 
we will show below, axinos are naturally produced at low 
temperatures of a few GeV, consistent with the Q-ball 
scenario of EMD but still before the period of Big Bang 
Nucleosynthesis (BBN). 

2. The Enqvist-McDonald scenario. 

We now briefly present the main features of the EMD 
variant of AD baryogenesis. It is assumed that the AD 
field (/) is a _D-flat direction in the MSSM. Its potential is, 
in general, lifted by soft supcrsymmetric (SUSY) break- 
ing terms and nonrenormalizable terms [17, 18]. 

The potential of the AD field, including inflaton- 
induced terms, reads 
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where is the soft SUSY breaking mass for the AD 
field and a radiative correction is given by ifln|(/)|^. A 
flat direction dependent constant, K, takes values from 
—0.01 to —0.1 [19, 20]. A denotes a renormalization scale 
and —ciH^, with ci ~ 1, is the negative mass-squared 
term induced by the energy density of the inflaton [18]. 
Terms proportional to A and C2 are the trilinear terms 
from low energy SUSY breaking and those induced by the 
inflaton, respectively, while TO3/2 denotes the gravitino 
mass. The nonrenormalizable terms in Eq. (1) come from 
the superpotential W = A/nM"~^(^", where A is the 
Yukawa coupling and M is some large scale acting as a 
cut-off. In SUGRA, it is natural to assume M = Mp ~ 
2.4 X 10^* GeV which is the reduced Planck mass. 

Since during inflation the Hubble parameter H ^ 
~ inj,i2, the AD field settles down at the minimum 
of the potential (1) which is given by 
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It is clear that the AD field can naturally develop a very 
large VEV, which is possible in nonminimal Kahler po- 
tentials [18], or if large enough trilinear term A is induced 
by the inflaton [21]. 

We have neglected in Eq. (1) thermal mass terms 
/i^T^ 1^1^, where h denotes couplings of the AD fleld to 
other particles [22]. They would play a role if the AD 
field VEV were relatively small. We have also neglected 
two loop thermal efi'ects due to the running of gauge 
coupling which generate a term aT^ln (|(^|^/T^), where 
|a| = 0(10^2) [23]. They will not be important below. 

As H decreases, the AD field traces the instantaneous 
minimum after inflation, begins to oscillate when H^^^ ~ 
rrv^ and, after a few dozen turns, produces a nonzero 
baryon number and then fragments into Q-balls. 



The baryon number density for the AD field (j) is given 
by Jib = iq{4>*(j) — 4>*(j)) where q is the baryonic charge for 
the AD field. By using the equation of motion of the AD 
field, the charge density can be rewritten as 
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with a{t) being the scale factor. When the AD field starts 
to oscillate around the origin, the baryon number density 
is induced by the relative phase between A and 02- With 
the entropy density after reheat s = 47r^g(*T^/90, we can 
express the baryon asymmetry as 
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Here, tosc denotes the time of the start of the oscillation, 
and sin 5 is the effective CP phase. 

From now on, we consider the case of n = 6 because a 
promising AD field for our scenario, a udd direction be- 
longs to this class. Let us first evaluate rib/s. The baryon 
asymmetry, Eq. (4), for the relevant case is estimated as 
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which is of the right order. (The previously made as- 
sumption M ^ Mp is crucial, for otherwise Q-balls would 
decay too early or would eveporate.) The low reheat 
temperature Tr after inflation is required to explain the 
appropriate baryon asymmetry. In this case, the AD con- 
densate fragments into Q-balls [19, 20]. 

The growth of perturbations of the AD field and its 
subsequent fragmentation into Q-balls crucially depends 
on the logarithmic correction to the (f)'^ mass term in 
V{(j)), Eq. (1). An essential requirement is that U(0) 
is flatter than quadratic, or that K <0 [20] . This can be 
achieved in SUGRA models with a nonminimal Kahler 
potential [19]. 

In order to discuss the evolution of Q-balls, first we 
briefly summarize their relevant properties in GRMSB 
models. The radius of a Q-ball, R, is estimated as ~ 
2/{\K\m'^) [19]. The charge is roughly given by Q ~ 

^wR^nbiti) ~ ^ttR^ (Hi/Hosc)'^ nb\toi,cJ where the suffix 
i represents the time when the spatial imhomogeneity 
becomes nonlinear, which can be evaluated as [24] 



Q ~ 6 X 10 



_3 2gl^lsin(5m3/2Mp' 



3/2 



A3/2 



m 



5/2 



1 X 10 



,20 



A3/2 



q\A\smS / m3/2 x /ITeVy^' 
' ' UOOGeV/ \ J ' 



3 



Unless Q > 0(10^*), Q-balls will evaporate before decay- 
ing [25]. For Q as in Eq. (6), Q-ball decay temperature 
is Trf ~ 1 GeV to 1 MeV [4, 26]. For example [26], 



Q 



1/2 



which is lower than the typical freeze out temperature of 
WIMPs, Tf ~ mx/24. Thus, the LSPs generated in Q- 
ball decays do not subsequently thermalize. Nor will the 
baryon asymmetry be washed out by sphaleron effects 
since T4 < T^w [4]. Note that Q-balls decay prior to 
BBN and thus do not spoil its successful predictions. 

In the EMD scheme, T/j also must be rather low. (This 
justifies neglecting thermal effects in Eq. (1).) In fact, 
unless Tij <^ 10^"^ GeV, Q-balls could thermaUze [4]. In 
order to preserve the fit, - ^dm relation, one needs to 
suppress the neutralino population from freeze-out. For 
this to happen it would be sufficient to assume Tr <^ Tf . 

It is easy to see why in the EMD scenario, the ratio 
flh/fl^ should be less than 1. The Q-ball is basically 
a huge "bag" of squarks. It decays predominantly via 
q ^ q + X- Thus, for one unit of a baryon number, at 
least > 3 units of nonbaryonic number density are 
created. (This number can be larger than 3 if one takes 
into account additional decays of squarks into heavier 
charginos and neutralinos which then cascade decay into 
the lightest neutralino, which arc model dependent.) In 
other words, the LSP number density after Q-ball 
decay is given by riy^^ = N-^^fsnt, where Jb is the frac- 
tion of baryon asymmetry carried by the AD field (p that 
is transferred into Q-balls. Prom lattice calculations, 
/b — 1 [27]. Assuming that the LSPs subsequently do 
not undergo any significant self-annihilation, and since in 
general flh"^ = mY = mn/s, this can be recast into 



rrir, 



Ofe _ TUnYb 



(7) 



where m„ denotes the mass of a nucleon and ni)^ the mass 
of the neutralino. It is clear that Eq. (7) implies fi^/fi^ to 
be less than one but not <C 1. In the EMD scenario, not 
only both types of matter are simultaneously produced 
but also a right ratio of their abundances is predicted. 

Unfortunately, this attractive picture runs into a se- 
rious problem of over-producing neutralinos, as noticed 
already by EMD themselves [5, 28]. Since Ym ~ 3.9 x 
10-1° (Q/iVO.ll) GeV, one can rewrite Eq. (7) as 
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In order to remain consistent with the values of Ub and 
fi^^/i^ derived from observations, the neutralino mass has 
to be 0(1 GeV) which, in the MSSM, is excluded by 
LEP [10]. Here, we have neglected a possible contribution 
to the LSP density from freeze-out. If it were significant, 
the problem would become only worse. Moreover, the 
condition (8) puts into question an attractive AD mech- 
anism in a large class of SUGRA models. 



To circiimvcnt these problems, one has to review as- 
sumptions in the above discussions, namely: 

(i) LSPs produced in Q-ball decay do not annihilate; 

(ii) The LSP is the lightest neutralino of the MSSM. 
If we relax assumption (i), the neutralino LSP with 

the mass of 0(10^ GeV), consistent with LEP, becomes 
acceptable. Indeed, allowing for subsequent LSP self- 
annihilation, the LSP density will be reduced by ~ 
{a^v)Td [29]. If the cross section cr^ for the LSP (sclf- 
) annihilation is large enough, e.g., when the LSP is Higgs- 
ino or Wino-like [29] , the relation between Qb and CIdm 
is lost. One interesting exception is when the energy den- 
sity of universe is dominated by Q-ball itself [30]. 

Alternatively, if we lift assumption (ii), the flb-^DM 
relation may be preserved. One way is to consider, e.g., 
models with the Higgs sector supplemented by a singlet. 
If its fermionic partner, the singlino, is the LSP then, 
for some specific choices of parameters [31], it could be 
possible to circumvent the LEP bound and perhaps also 
to suppress the LSP abundance from freeze-out. In the 
rest of the paper, we will investigate axino LSP as DM. 

3. Axino dark matter from Q-balls. 

In general, axinos, like gravitinos, can be produced in 
both thermal processes (TP) and in nonthermal processes 
(NTP), e.g., in late decays. TP consists of the scatterings 
and the decays of particles in the thermal bath. NTP 
is given by the decay of the Next-to-LSP (NLSP) relic 
(which, for simplicity, we assume to be the neutralino) 
from freeze-out or from the decay of Q-balls in our sce- 
nario. 

The relevant Boltzmann equations can be written as 

% + SHn^ = -{a^v){nl - nl^^^) + Jq - T^n^p) 
fia + SHria = {<Tv)ijninj + {uv)ini + Ty^n-y^, (10) 

where (cp^w) is the usual neutralino freeze-out term, 7q 
denotes the contribution to x from Q-balls decay, P^^ is 
the decay rate of the neutralino, {av{i + j a + 
and {av{i ^ d+ are the scattering cross section and 
the decay rate for the thermal production of ax;inos. 
The total NLSP abundance is given by 



NxfB^+Y' 



(11) 



where Y^^ ~ H/si^^^^m^/Tf/ {av)ann, as usual. Since 
ns, — n-x, owing to R-parity conservation, the resulting 



number density of axino is given by 

V — vNTP , T^TP 
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(13) 



where Y^^ denotes the axinos produced by thermal pro- 
cesses. Since typically Y^^ ^ lO-n, one can sec that 
nonthermal production of axinos due to the thermal relic 



4 



NLSPs decay can easily be negligible compared to that 
from Q-ball production, and its contribution to i^a/i^ 
is further suppressed by the ratio ma/m^. The ther- 
mally produced axino Y^^ also can be subdominant, say 
< 10-", for Tr ^ 100 GeV [15, 32]. Hence, the ax- 
ino dark matter density is estimated as 




One can see that the baryon asymmetry and the dark 
matter abundance are readily linked. 

As mentioned above, axino mass is strongly model de- 
pendent; in particular, it critically depends on how the 
visible and hidden sectors are coupled [11, 12, 13]. At tree 
level, either ma = 0(777.3/2) or 0{m'^^^/ fa) = O(keV). 
However, in the latter case, trilinear terms can generate 
a substantial 1-loop correction of order Jq/Stt'^A, where 
/q is the Yukawa coupling of the heavy quark to a singlet 
field containing the axion, which gives ma in the range 
of a few tens of GeV or less [12, 13]. 

The final check point is the compatibility with suc- 
cessful predictions of BEN. However, this is not really 
a problem for the axino LSP because its interactions 



are less suppressed than those of the gravitino, roughly 
by {Mp/fa f and, so long as the NLSP is heavier than 
about 150 GeV, axinos are produced before the time of 
BBN [15]. In contrast, the gravitino LSP would be pro- 
duced in late NLSP neutralino decays, which faces strong 
constraints from BBN [33, 34]. 

4. Conclusions. 

We have shown that the framework with cold dark 
matter axino LSP produced in Q-ball decays can explain 
the abundance of dark matter and the baryon asymmetry 
simultaneously and may be an answer to the fib ~ floM 
coincidence. In this approach, the similarity between Qb 
and ft DM is explained by basically only the axino mass of 
order of (sub-)GeV. The essential property of Q-ball de- 
cays is that one can predict the number of SUSY particles 
per one baryonic charge from Q-ball A^^^. A characteristic 
feature is low reheat temperature Tr of 10^ GeV. 
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